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i-O Abstract. The constraint that can be set on the Hot 1. Introduction 
Dark Matter (HDM) model of structure formation from 
the temperature distribution function (TDF) of X-ray 
clusters of galaxies is examined in terms of the ampli- 



tude to be given to the primeval matter fluctuations for 
consistency with present observations. 

The TDF of clusters is derived for different normaliza- 
I tions of the HDM power spectrum in an Einstein-de Sitter 
£SJ universe (fio = 1) using the statistics of peaks of a random 
Gaussian field. 

Concerning cluster formation only, the neutrino pic- 
ture fails to reproduce the TDF for Ho = 50 km s _1 
45 Mpc -1 but appears to be marginally consistent with ob- 
Onservations for Ho = 70 km s _1 Mpc -1 provided the nor- 
Q malization is a « 0.42, in terms of the ratio a = <tq/6 s 
with (To the rms variance of the fluctuations linearly 
(Z2 evolved to the present epoch and 6 S the collapse threshold 
^ required in structure formation models. 

However, using the constraint that the formation 
epoch of QSO's puts on the HDM model, it appears 
that this normalization is too small to be consistent with 
quasar formation. Allowing higher values for the Hubble 
constant can alleviate this problem but leads to an un- 
acceptable short age of the universe. Moreover, the con- 
straint put by the recent measurement of the level of fluc- 
tuations in the cosmic background radiation by the COBE 
DMR experiment (Smoot et al. 1992) tends to reject mod- 
els with high values of Ho- Despite the fact that the ob- 
servational constraints do not appear as stringent as they 
were claimed to be in some previous analysis, we still con- 
clude that the neutrino model fails in any case to account 
for the whole set of available observations in a consistent 
way. 
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Much work has already been done on the once popular 
Hot Dark Matter (HDM) model of structure formation 
(e.g: Cowsik k McClelland 1973; Bond et al. 1980; Peebles 
1982; Bond k Szalay 1983; Kaiser 1983; Frenk et al. 1983; 
White et al. 1983, 1984; Melott 1985, 1987; Buchert 1991; 
Zeng k White 1991). Several attractive features naturally 
arise in this scenario: 

1) first of all, while whether they are massive or not is still 
an open question, neutrinos do exist. Stringent limits on 
the electron neutrino mass have been obtained by differ- 
ent methods but the constraints put on the muon and tau 
neutrino masses do not exclude a mass of a few tens of 
eV. While the standard Big Bang nucleosynthesis models 
predict that the total amount of baryonic material in the 
universe should not exceed about one tenth of the closure 
density of the universe, the existence of a non baryonic 
dark matter constituent like one species of massive neu- 
trinos might reconcile the abundance of baryonic material 
predicted by the nucleosynthesis with the exigence of a 
density parameter fio close to 1 as expected in inflation- 
ary scenarios. 

2) the short-wavelength cut-off in the HDM power spec- 
trum provides a natural explanation to the filamentary as- 
pect of structure on large scale (de Lapparent et al. 1986) 
and to the large coherence observed in the real galaxy 
distribution. 

However, the model seems not to resist to a closer in- 
vestigation. As pointed out by Peebles (1982), then em- 
phasized by the numerical experiments of White et al. 
(1983, 1984), the neutrino picture suffers from a timing 
problem that comes to light when comparing the large- 
scale structure properties at the present epoch with the 
small-scale constraint at high redshift. This problem can 
be conveniently expressed in terms of the normalization 
of the model through (To, the rms variance of the matter 
fluctuations linearly evolved until now. Indeed, each ob- 
servational constraint can be related to a value required 
for (To. In general, a single observational constraint will 
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be satisfied for a certain range of values of (To. A robust 
test of the model is obtained by examining simultaneously 
several constraints, i.e. the consistency of the values of (To 
required by the various observations. In their pioneering 
work, White et al. (1984) choose the onset of galaxy for- 
mation as the time when one percent of the mass in their 
simulations has undergone collapse. This happens, in their 
simulations, after an expansion factor of 3, when the linear 
rms fluctuation is a = 0.66. In order to take into account 
the existence of high redshift quasars (z > 2), they nor- 
malize their simulations to (To = 0.66(1 + z c ) ps 2.2 and 
(Jo ~ 4.4 so that the first collapse of objects occurs at 
z c = 2.5 and z c = 6 respectively. As a result of using such 
high normalizations, the clusters produced later in the 
simulations, i.e. at z = 0, have a much higher mass (and 
temperature) than the one observed in the real galaxy 
distribution. From the argument that HDM cannot repro- 
duce consistently both the small-scale, i.e QSO's, and the 
large-scale, i.e. clusters, properties of structure, the model 
was early ruled out (White et al. 1984; Kaiser 1983). How- 
ever, it should be emphasized that the choice of the onset 
of QSO's formation is quite uncertain as they represent a 
negligible fraction of the total collapsed material and are 
below the resolution limit of White et al. numerical exper- 
iments. In retrospect, their criterion may seem too drastic 
and their conclusions have thus been questioned (Melott 
1985; Melott 1987; Centrella et al. 1988); accordingly, the 
model has been revisited from time to time (Anninos et 
al. 1991; Zeng k White 1991; Cen k Ostriker 1992; Tuluie 
et al. 1994). 

Recently, Blanchard et al. (1993) have revisited the 
constraints set from small-scale (QSO's) observations in 
the neutrino picture. Their semi-analytical approach leads 
to the conclusion that, considering the large uncertainty 
that exists on the density threshold 6 S involved in struc- 
ture formation models, the existence of high redshift 
quasars can be taken into account if the normalization 
of the spectrum, (To, ranges from 1 to about 3. The nor- 
malizations used in the pioneering work of White et al. 
(1984) clearly do not cover all this range. 

At this point, it seems interesting to us to revise 
the constraint that can be set on the HDM model from 
galaxy clusters. In this paper, we propose an analytic ap- 
proach using the cluster temperature distribution function 
(TDF), i.e. the abundance of galaxy clusters as a function 
of their X-ray temperature, to constrain the normaliza- 
tion of the HDM spectrum. The constraints that can be 
set on galaxy formation theories from the TDF of clusters 
have been discussed in various models (Henry k Arnaud 
1991; Blanchard k Silk 1991; Bartlett k Silk 1993; Ouk- 
bir k Blanchard 1992), while some authors have rather 
used the constraint imposed by the velocity distribution 
function (Evrard 1989, Lilje 1990). However, this latter 
method suffers from uncertainties as cluster velocity dis- 
persion might be overestimated (Frenk et al. 1990) be- 
cause of projection effects along the line of sight. More- 



over, the observed galaxy velocity dispersion, a ga i, might 
not trace the matter velocity dispersion: there might ex- 
ist a velocity bias as claimed by Carlberg k Couchman 
(1989) according to their simulation results. On the con- 
trary, X-ray data appear to be more reliable (Edge et 
al. 1990). The X-ray luminosity function has been also 
used to constrain galaxy formation models (Kaiser 1986, 
Evrard k Davis 1988). However, the luminosity is entirely 
dominated by the core properties of the X-ray gas that has 
been shown to reside in clusters and it is unlikely that its 
scaling properties can be modeled confidently (Blanchard 
et al. 1992b). On the contrary, the gas temperature is a 
quite direct measure of the potential depth and the X-ray 
temperature distribution function appears as a clean test. 

In our approach, the main hypothesis is that the clus- 
ter gas is close to isothermal hydrostatic equilibrium at 
least out to the virial radius. Observations tend to con- 
firm isothermality in the central part (apart from clusters 
in which evidence of cooling flows has been found). Hydro- 
dynamical simulations from Evrard k Davis (1988) tend 
also to favour this hypothesis and we adopt it hereafter. 

Our derivation of the cluster TDF is based on the 
statistics of peaks of the initial matter density field (see 
Bardeen et al. 1986, BBKS throughout this paper) and 
on the spherical collapse model (Lemaitre 1933; Gunn 
k Gott 1972; Peebles 1980). This method appears well 
designed for an HDM scenario since peaks of the field 
expected to lead to structure formation are unambigu- 
ously defined (because of the absence of structure on small 
scales), whilst the same method for the Cold Dark Mat- 
ter model, as for any hierarchical scenario, faces the so- 
called 'cloud-in-cloud' problem when the field is smoothed 
from small scales to larger scales (BBKS, Peacock k Heav- 
ens 1985). Using the spherical infall model is however 
more questionable as in a 'pancake' process the collapse 
of matter is expected to proceed along one privileged axis 
leading to the formation of highly anisotropic structures. 
The collapse of peaks might occur for a density thresh- 
old 6 S smaller than the canonical value 1.68 derived in 
the top-hat model. However, high-j^ peaks are expected to 
be spherical (BBKS). Moreover, it has been shown that 
high-j^ density fluctuations rigorously follow the spherical 
collapse for any cosmological model (Bernardeau 1994), so 
that this value, <5 S = 1.68, seems to be fair. Nevertheless, a 
possible way to get rid of the systematic uncertainty on 
(To induced by the uncertainty on the threshold 6 S is to ex- 
press the derived normalizations in terms of the parameter 
a (the inverse of the lowest v of structures that form): 

a = <r /6 s 

rather than merely <tq. Blanchard et al. (1993) pointed out 
that the quantity a is directly constrained by the obser- 
vations independently of the uncertainty on 6 S . This will 
allow us to compare in a consistent way the normalization 
we obtain on cluster scale to the one derived by Blanchard 
et al. for QSO's. 
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2. Cosmological framework and fluctuation spec- 
trum 

We consider the growth of large scale inhomogeneities of 
the initial matter density field in an Einstein-de Sitter 
(flo = 1) universe according to gravitational instability 
theory (in the following, the subscript for the density 
parameter and the Hubble constant refers to present val- 
ues). The statistical properties of the field, assumed to be 
random and Gaussian, are described by P(k, t), the power 
spectrum of the fluctuations, i.e., the correlation function 
of the Fourier transform of the density contrast 6 defined 
as the local excess of matter over the mean density of the 
universe: 



P(M) = (\6(k,t 



6(r) 



P( r ) ~ (p) 
(P) ' 



(1) 



In the HDM model, the assumption is made that the 
universe consists mainly of massive neutrinos and in an 
fio = 1 cosmology, according to primordial nucleosynthe- 
sis, neutrinos must then account for 95% of the total mat- 
ter if Ho=50 km s _1 Mpc -1 . Such an assumption con- 
strains the power spectrum except for its normalization 
(the amplitude of the primordial fluctuations) which is 
the only main free parameter of the theory. It is usual to 
choose the Harrison-Zel'dovitch scale invariant spectrum 
as the initial spectrum: P,(k, i) oc k n with n = 1. The 
power spectrum at recombination is obtained by multiply- 
ing Pi by a transfer function that accounts for the history 
of the fluctuations (free streaming damping in the neutrino 
hypothesis). We choose the transfer function derived from 
Bond & Szalay (1983) and given by BBKS. The inferred 
spectrum is characterized by a cut-off A c which depends 
only on the density parameter of the neutrino constituent 
fi;, (for one species) and the Hubble constant Ho = lOO/i 
km s -1 Mpc -1 : 



A, 



13 (/i 2 ^) -1 Mpc. 



We approximate S]„ « 1. The spectrum presents a lack 
of power at small scales, which prevents small structures 
such as galaxies from forming first. Clusters then appear 
as the first structures to form in this model. 

3. Temperature distribution function of galaxy 
clusters 

In the gravitational instability picture the maxima of the 
matter density field are likely to be the privileged sites 
of collapse of non-linear objects such as clusters of galax- 
ies (BBKS 1986). In hierarchical models, the counting of 
peaks is not reliable because of the 'cloud-in-cloud' prob- 
lem and a Press & Schechter-like approach turns out to be 
more adapted (Blanchard et al. 1992a). In the HDM sce- 
nario, contrary to hierarchical scenarios, there is one def- 
inite epoch at which the first structures form. This epoch 
corresponds to the collapse of the highest peaks of the 



field. It is only much later, when (To ^ 1, that fluctuations 
on large scales will collapse, embedding the first generation 
objects. In this scenario, the count of non-linear objects 
from the peaks of the density field is thus a reliable ap- 
proach at least during the collapse of the first objects (see 
below for a more quantitative discussion of the 'cloud-in- 
cloud' in the HDM model). During the collapse of the first 
peaks, the gas is shock-heated, reaches the virial temper- 
ature and is gravitationally confined through an hydro- 
static equilibrium. The final temperature depends on the 
total collapsed mass. In order to identify clusters hotter 
than T, i.e. more massive than some value M, we first 
smooth the density field by a top-hat window function of 
radius R (with M = p^irR 3 ) and we assume that clusters 
form at location where the smoothed field is above the 
threshold 6 S . The choice of an adequate filtering window 
for smoothing the field is one difficulty in this approach. 
The spherical (top-hat) window is probably the most rele- 
vant if the collapse of structure is expected to be spherical. 
However, the sharp features of the window in real space 
induces oscillations in Fourier space, which adds power to 
the spectrum out to several cut-off scales. The Gaussian 
window avoids such a difficulty but the relation between 
the mass of a collapsed object and the mass of the cor- 
responding peak in the field is not obvious in this case 
(Peacock & Heavens 1985). However this problem is of 
minor importance for HDM since the low-pass filtering 
procedure induces high-frequency spikes on scales where 
there is no power. We therefore expect that Gaussian and 
top-hat filtering lead to comparable results in this case. 

The temperature distribution function Af(T) of clus- 
ters at some epoch z v can then be directly related to the 
statistics of peaks of the smoothed field: 



Af(> T) 



Af pk {v, R k )dv 



M(f)df (2) 



In equation (2) we use the formalism developed by BBKS, 
so that Afpk, v and are defined by: 



6 



:<7i 



v= — (l + z„), R* = a/3- 



(3.2) 



where G(y,yv) in equation (3.1) is the function given by 
equation (A. 19) of BBKS and the parameter y is defined 
by: 

(3.3) 



C2C0 



7 and R* are called the spectral elements of the field and 
they involve the spectral moments a, defined below: 



-fit) 



k 2 dk 
2tt 2 



P(k,t)k 2i . 



(4) 



ao(Rf) is the rms fluctuation of the field filtered at the 
comoving scale Rf. ai(Rf) and <T2{Rf) are related to the 
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first and second derivative of the field. Mpk{v) is the mean 
differential number of peaks per comoving volume unit (cf 
BBKS). can be related to the mean distance between 
peaks in the following way: the cumulative number density 
of peaks of arbitrary height is: 



ipk(-oo) 



M pk (v)dv ps 0.016i?" 3 . (BBKS : 4.116) 



The mean distance between peaks of arbitrary height is 
then: 

< d ph >= (0.016i?^ 3 )"^ ~4x^. 

The approach considered in equation (2) has been 
known to be plagued, in a hierarchical scenario of struc- 
ture formation, by the so-called cloud-in-cloud problem 
(BBKS 1986, Peacock & Heavens 1985). This term re- 
covers the fact that two (or more) separate peaks in the 
field smoothed at some scale R\ might be part of a single 
peak, i.e. structure, at some other scale R2 > R\. Such an 
approach is therefore likely to overestimate the real num- 
ber of structures in a hierarchical scenario since several 
non-linear peaks fated to end up into a single structure 
might be counted as individual objects. This is probably 
the major issue that makes the 'peak approach' a crude 
way of estimating the temperature function in this kind 
of scenario. 

In the case of a top-down type model as the HDM 
model envisaged here, because of the lack of power on 
small scales, the mean distance between peaks is large 
enough at all the filtering scales of interest, i.e. on cluster 
scales, to avoid the possibility that separate peaks on one 
scale are absorbed in a unique structure on a larger scale. 
The mean distance between peaks of arbitrary height is 
4 x 13.69 ~ 55 Mpc (h = 0.5) for HDM and the mean 
distance between peaks on cluster scale (8/i _1 Mpc) is 
~ 93 Mpc. We thus expect that peaks in the HDM scenario 
are unambiguously defined and that equation (2) gives a 
good approximation to the number of non-linear objects. 

Deriving this equation leads to: 



M(T) = M pk (v s ,R+)'^- 



dM ph {y,R*) dR* 

dR* dl/] dr" 



(5) 



The first term in the second member of equation (5) counts 
the number of peaks with height v s within dv s while the 
second term counts the variation in the number of peaks 
above the threshold 6 S when the smoothing scale changes 
so that R* increases of dR*. As it will be clear on figure 1 
below, the second term can be neglected on cluster scales 
so that equation (5) can be approximated by: 



dv s 



which can be rewritten, according to (3.2), as: 



— 1 dao 
rf~dT' 



(6) 



(7) 



Note that this equation depends on a = <j /S s and not on 
(Jo or 6 S separately so that it is not necessary to specify the 
value of 6 S at this level. In equation (7), <tq is related to 
the temperature through the smoothing scale Rth and: 



dcrp 
~dT 



dcrp 
dR 



dR 



TH 



TH 



dT 



(8) 



The scale Rth is related to the mass of the peak by: 



M 



4ir . 



-R TH p 



(9) 



where p is the mean density of the universe and a rela- 
tion between the mass and the temperature of a newly 
formed object (i.e. galaxy cluster) is needed to evaluate 
dRxH /dT. We derive this relation from the assumption 
that the intracluster medium, the baryonic gas that has 
been shown to reside in clusters, is in an hydrostatic equi- 
librium state at constant temperature. Such an assump- 
tion is discussed below. In this case, the equation of hy- 
drostatic equilibrium for the gas is merely: 



k R T. 



—pm p GM 
Jp + JT R 



(10) 



dlnpg 
dlnR ' 



It 



dlnT 
dlnR 



where p is the mean molecular weight of a fully ionized 
plasma (p ps 0.63), ks is the Boltzmann constant, m p is 
the proton mass, p g is the intracluster gas density and 
M is the mass enclosed in a spherical shell of radius R. 
Equation (10) leads to: T oc M 1 2 5 /3 (l + z v )h 2 ' 3 keV with 
z v the epoch of formation and Mi 5 the cluster spherical 
mass, related to the smoothing scale as mentioned above 
(equation (9)), in unit of 10 15 solar mass. We normalized 
this relation to Evrard's hydrodynamic simulations result 
(Evrard 1990) and we obtained: 



T ps 6.4M 1 2 5 /3 (1 + z v )h 2/3 keV. 



(11) 



Evrard's numerical simulations have been performed 
for a CDM initial power spectrum and using them for the 
neutrino picture is questionable. However, the derivation 
of (11) involves mainly the equation of hydrostatic equilib- 
rium and the expression of the mass enclosed in a sphere 
of A p» 200, the density contrast of a virializing spheri- 
cal cluster according to the spherical collapse model. As 
these two relations are independent on the characteristics 
of the underlying fluctuation spectrum, and as cluster col- 
lapse appears to be insensitive to the amount of small-scale 
power present in the spectrum (see Evrard & Crone (1992) 
simulations), we expect (11) to remain valid in the neu- 
trino scenario. Lilje (1990) made the same assumption. Let 
us mention again that equation (11) assumes that the gas 
in galaxy clusters can be regarded as isothermal. This as- 
sumption is favoured by observations (Hughes et al. 1988; 
Hughes 1992) but is uncertain far from cluster cores. 
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The only remaining quantity to be determined in equa- 
tion (7) is the ratio of the overall amplitude of the fluctua- 
tions ((To at Mpc), to the threshold 6 S , i.e. the parameter 
a. Note that (To, which is nothing but an integral over the 
fluctuation spectrum according to definition (4), diverges 
in for CDM but converges for HDM. We propose to 
determine the quantity a by scaling relation (7) on the 
observed temperature function of galaxy clusters given by 
Edge et al. (1990). This is one of the main purpose of this 
paper. 

4. Results and discussion 

We present on figure 1 the temperature distribution func- 
tion of galaxy clusters as obtained from equation (5) (solid 
line) and from equation (6) (dashed line). The curves are 
plotted for different redshifts of cluster formation. The 
normalization of the spectrum is arbitrary and has been 
taken to be a 0.36. The graphs are given for h = 0.5 
which is required in order to obtain an acceptable age 
of the universe. Because of a normalization in the linear 
regime, strong evolution effects between redshift 1 and 
now can be noticed. The main feature of the temperature 
distribution function is a plateau for temperatures below 
1 keV, which corresponds to galaxies or small groups. This 
feature clearly results from the lack of power at small 
scales in the neutrino picture spectrum. In Figure 2, we 
present the temperature distribution function for different 
normalizations of the spectrum at z v = obtained from 
equation (6), compared to observations. The filled squares 
correspond to Exosat and Einstein observatory data on a 
sample of 55 X-ray clusters (Edge et al. 1990). The filled 
square at 13.6 keV is from M. Arnaud et al. (1992). Er- 
ror bars are given at 90% confidence level. In order to 
fit correctly the cluster abundances with temperature in 
the range 5 to 15 keV, it is necessary to have the am- 
plitude a in the range 0.33 to 0.39 (these bounds do not 
depend on the actual value of the threshold of collapse 6 S ). 
Higher amplitudes would copiously overproduce high tem- 
perature clusters. On smaller temperatures (T < 5 keV), 
the curve conflicts with data (points PI and P2 on figure 
2a). The observed number density of clusters at T 2.2 
keV (point PI) is k, 2.7 10" 6 Mpc" 3 while the total num- 
ber density of clusters we derive at this temperature from 
the peak approach is « 1.5 10 -7 Mpc -3 . In P2 (T « 3.5 
keV), the observed number density is k, 1.4 10 -6 Mpc -3 
while we find k, 1.4 10" 7 Mpc" 3 . For HDM, the total 
number of peaks (i.e. the number of clusters as in this sce- 
nario one peak leads to the formation of one cluster) with 
height > v = a -1 k, (0.35) -1 k, 2.85 is (see figure 3 from 
BBKS with 7 k, 0.74 and k, 13.69 Mpc): k, 6.2 10" 7 
Mpc -3 . The total observed number density of clusters for 
PI and P2 is k, 4.1 10" 6 Mpc" 3 , about one order of mag- 
nitude higher than the total number density of clusters 
expected in the HDM scenario so that the HDM model 
can not reproduce the observed TDF for h = 0.5. How- 



h=0.5 _ 




1 10 
T [keV] 



1000 



Fig. 1. Temperature distribution function from the peak ap- 
proach for different redshifts. The spectrum is arbitrarily nor- 
malized to a fa 0.36. The solid line is obtained from equation 
(5) and the dashed line is from equation (6). The second term 
in equation (5) turns out to be negligible. As the spectrum 
is normalized in the linear regime, strong evolution effects are 
noticed between redshift 1 and the present epoch. The main 
striking feature of the TDF is a plateau for temperatures below 
1 keV. 



ever, temperatures lower than 5 keV correspond to groups 
or small clusters scales which might form by fragmenta- 
tion of pancakes in the HDM model. Structures formed by 
fragmentation of peaks have not been taken into account 
in our analysis and fragmentation processes may then ex- 
plain the discrepancy between the curve and the data. We 
can test this possibility by the following procedure. Frag- 
mentation may increase the number of clusters and as the 
'new' clusters will be in the neighborhood of the 'parent 
peak', this will translate into a high correlation between 
clusters. From the observations, we can evaluate the mean 
number of clusters N v associated with a collapsed struc- 
ture, i.e. the mean number of neighbours of a given cluster: 



N„ 



(12) 



where n c is the comoving number density of clusters, S, c {r) 
the cluster-cluster two-point correlation function, and R 
the comoving radius of the associated peak in the field. 
The cluster-cluster correlation function reads: 



(13) 




Fig. 2. a-c Temperature distribution function from the peak approach for different normalizations of the HDM spectrum 
compared to Edge's (1990) X-ray data: a h = 0.5, b h = 0.7, c h = 1. The cluster at 13.6 keV is from M. Arnaud et al. (1992). 
Error bars are given at 90% confidence level. The arrows on 2. a denote the upper limits to which the solid curve can be pushed 
up in the presence of fragmentation. 



with r pa 25/} -1 Mpc and j Pa 1.8 (H auser & Peebles 
1973). Using (13), equation (12) leads to: 

n, = ^ - n \ _ 7 R 3 r f (14) 

(3 - 7 )r 7 

where the radius R has been chosen to be roughly R+, 
which tends to overestimate N v . For point PI (figure 2a), 
n c is estimated to be 2.2 x 1.2 10" 6 Pa 2.7 10" 6 Mpc" 3 
at T « 2.2 keV and (14) leads to N v pa 1.1 for h = 0.5. 
For P2, we estimate n c pa 1.4 10 -6 Mpc -3 and we find 
N v pa 0.67. In the presence of fragmentation, the actual 
number of peaks is: 

M pk x (1 + N V ) Mpc" 3 . (15) 

On figure 2a we put the upper limits to which the solid 
curve can be pushed up due to fragmentation processes. 
Fragmentation turns out not to be efficient enough to ex- 
plain the discrepancy between the model and the obser- 
vations. Note that the low-end slope of our model TDF 
should be greatly in error to reconcile the curve and the 
data. 

We therefore conclude that the HDM model with 
h = 0.5 fails to reproduce the TDF of clusters. On fig- 
ure 2b and 2c, the same graphs as figure 2a are presented 
for h = 0.7 and h = 1 respectively. While the normaliza- 
tion appears to be only slightly sensitive to the Hubble 
parameter, the fit is improved in the case of h = 0.7. For 



h = 1, the fit is no longer good and the age of the uni- 
verse is unacceptably small. We find that for a pa 0.42 
HDM with h = 0.7 reproduces roughly the observed clus- 
ter temperature function. 

Figure 2 shows also that the TDF of clusters is very 
sensitive to the initial amplitude of matter fluctuations: 
increasing the normalization by a factor 2 results in about 
a factor 4 in temperature. In other words, the comparison 
between the TDF and observations provides a powerful 
constraint on a. We draw on figure 2a and 2b the TDF 
for a = 0.59. For 6 S = 1.68 this corresponds to a normal- 
ization (To Pa 1. It appears clearly that this normalization 
is inconsistent with the data (for 6 S < 1.68, the discrep- 
ancy is reinforced). A fortiori, <tq values used by White 
et al. (1984) (cf introduction) are, as expected, much too 
large to account for cluster formation consistently with 
observations. Note that, as a check of our calculation, the 
number density of neutrino clusters we obtain at 40 keV 
with their normalizations is in good agreement with their 
simulation results. 

The Press & Schechter (1974) formalism is not ex- 
pected to provide a good estimation of the mass func- 
tion in the HDM model. Specifically, Press & Schechter 
assumed that the total mass in non-linear objects with 
mass M is related to the probability of finding a sphere 
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that satisfies the non-linear criterion: 

1 



P(> v.) 



'2tt 



-e 2 dv, 



(16) 



the field having been smoothed on the scale Rq = 
(3M/4tt)-. This measures only the volume correspond- 
ing to the centers of such non-linear spheres while the full 
volume of spheres of radius R > Rq should be taken into 
account; therefore (16) is an underestimate. The factor 2 
introduced by Press & Schechter (1974) to correct this ef- 
fect is purely phenomenological. However, as the Press & 
Schechter mass function is known to fit rather well the 
results of numerical simulations, we have checked to what 
normalization it leads. Using equation (11) and the Press 
& Schechter prescription for the mass function, we get: 



P M 1 + Z v) -(v 2 l2) 



7T M 



d<j 



dT 



(17) 



This approach leads quantitatively to the same results as 
those on figure 2. We find, in order to fit Edge's data: 
a pa 0.47 (h = 0.7). This result is slightly higher than the 
normalization derived previously. This may reflect the fact 
that the Press & Schechter approach is likely to underes- 
timate the true number of non-linear objects on the high 
mass end. In both cases however, the normalization ob- 
tained shows that according to HDM, most of the matter 
in the universe is still in the linear regime by now. 

Until now, we have shown that the neutrino picture 
is able to reproduce consistently with observations the 
galaxy cluster distribution in an universe with fio = 1 and 
h = 0.7 provided the normalization is adequately chosen. 
The validity of the model can be checked, as mentioned in 
the introduction of this paper, by comparing this normal- 
ization to the one obtained from QSO's by Blanchard et 
al. (1993). Recall that comparing directly the amplitude 
of a for clusters and QSO's cancels the uncertainty due 
to the unknown threshold 6 S . It appears from the con- 
straints put on quasar formation that «qso's has to be 
greater than 1.7 (h = 0.5). In this case, the ratio of the 
normalizations in terms of a is: 



«QSO's 
^Clusters 



4.8. 



(18) 



For higher values of the Hubble constant, this discrepancy 
can be reduced. For h = 0.7 and h = 1, «qso's ~ 1-4 and 
1 respectively so that: 



«QSO's 
^Clusters 



3.3, 



«QSO's 
^Clusters 



1.7. 



(19) 



For h = 1, the discrepancy is much smaller but the two 
values are still inconsistent. However if we use the most 
conservative assumption on 6 S for QSO's, i.e. 6 S 1 while 
clusters follow the spherical infall model (6 S = 1.68), the 
normalizations in terms of <tq become consistent for h = 1. 




Fig. 3. The different constraints put on the HDM model plot- 
ted in the a — h space. The cluster constraint is derived from 
the present work while the QSO's constraint is obtained from 
Blanchard et al. (1993). The COBE constraint is evaluated 
from Gorski et al. (1994) and Holtzman (1989). The hashed 
areas represent the allowed regions in the a — h space accord- 
ing to the corresponding constraint. 



But the age of the universe in this case would be only ap- 
proximately 6.5 Gyrs. 

We summarize our results on figure 3 where the dif- 
ferent constraints on the HDM model have been plotted 
in the a-h space. The cluster constraint has been de- 
rived from the work presented in this paper while the 
quasar constraint has been obtained from Blanchard et al. 
(1993). The COBE constraint is evaluated from Holtzman 
(1989) using the rms quadrupole coefficient value (Qrms 1 
k, 20 ± 5 yuK(2cr)) obtained by Gorski et al. (1994) for a 
spectral index n = 1. For the cluster constraint, we allow 
the virial temperature to differ from the observed value 
by a factor of two so that the allowed range for the nor- 
malization a is substantially increased with respect to the 
range obtained previously. Even with this wide range of 
allowed values for a, figure 3 shows that the HDM model 
in an Einstein-de Sitter universe can not account simulta- 
neously for the different observational constraints. 



Qr 



d2 where ci2 is the quadrupole coefficient in 



the spherical harmonic expansion of the sky temperature 
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5. Conclusion 

Concerning cluster formation only, we draw the following 
conclusions on the neutrino model: 

1) the HDM model in an Einstein-de Sitter (flo = 1) uni- 
verse with ft = 0.5 fails to reproduce the TDF. 

2) HDM in an fio = 1 universe with ft = 0.7 is marginally 
consistent with the available observations provided that 
the normalization is adequately chosen: a pa 0.42 (<tq pa 
0.7 for 6 S = 1.68). (Note that the age of the universe in 
this case is also marginally acceptable). 

3) HDM in an Einstein-de Sitter universe with ft = 1 is 
excluded by the temperature function constraint. (The age 
of the universe is also unrealistically low). 

Concerning small-scale constraint, Blanchard et al. 
(1993) concluded that QSO's can be accommodated with 
a lower normalization than those used by White et al. 
(1984). However, we find that the timing problem still re- 
mains as both large scale and small scale constraints can 
not be satisfied simultaneously for a single amplitude of 
the primordial fluctuations. Higher Hubble constants re- 
duce the discrepancy but lead to a very short age of the 
universe. Moreover, the level of CBR fluctuations as mea- 
sured recently by COBE tends to exclude a model with a 
high Hubble constant. We conclude that in any case the 
neutrino picture has to be ruled out. 

Acknowledgements. We thank the unknown referee for con- 
structive comments that helped us to improve the manuscript 
of this paper. We also thank Jim Bartlett for helpful discus- 
sions. C.B. acknowledges finantial support from the French 
Ministere de la Recherche et de l'Enseignement Superieur. 

References 

Anninos P., Matzner R.A., Tuluie R., Centrella J.M. 1991, ApJ 
382, 71 

Arnaud M., Hugues J. P., Forman W., Jones C, Lachieze-Rey 
M., Yamashita K., Hatsukade I. 1992, ApJ 390, 345 

Bardeen J.M., Bond J.R., Kaiser N., Szalay A.S. 1986, ApJ 
304, 15 

Bernardeau F. 1994, ApJ 427, 51 
Bartlett J.G., Silk J. 1993, ApJ 407, L45 

Blanchard A., Silk J. 1991, Proc. 26th Rencontres de Moriond 
in Astrophysics (Gif-sur-Yvette: Editions Frontieres) p. 93 

Blanchard A., Valls-Gabaud D., Mamon G.A. 1992a, AfeA 264, 
365 

Blanchard A., Wachter K., Evrard A.E., Silk J. 1992b, ApJ 
391, 1 

Blanchard A., Buchert T., Klaffl R. 1993, AfeA 267, 1 
Bond J.R., Efstathiou G., Silk J. 1980, Phys. Rev. Lett. 45, 
1980 

Bond J.R., Szalay A.S. 1983, ApJ 274, 443 

Buchert T. 1991, In: Physical Cosmology, eds: Blanchard A. et 

al. , Frontieres Paris, pp. 475-483 
Carlberg R.G., Couchman H.M.P. 1989, ApJ 340, 47 
Cen R., Ostriker J. P. 1992, ApJ 399, 331 

Centrella J.M., Gallagher III J.S., Melott A.L., Bushouse H.A. 
1988, ApJ 333, 24 



Cowsik R., McClelland J. 1973, ApJ 180, 7 

Edge A.C., Steward G.C., Fabian A.C., Arnaud K.A. 1990, 

MNRAS 245, 559 
Evrard A.E., Davis M. 1988, Nature 333, 335 
Evrard A.E. 1989, ApJ 341, L71 
Evrard A.E. 1990, ApJ 363, 349 
Evrard A.E., Crone M.M 1992, ApJ 394, LI 
Frenk C.S., White S.D.M., Davis M. 1983, ApJ 271, 417 
Frenk C.S., White S.D.M., Efstathiou G., Davis M. 1990, ApJ 

351, 10 

Gorski K.M., Hinshaw G., Banday A.J., et al. , 1994, Sissa 

Preprint, astro-ph/9403067 
Gunn J.E., Gott J.R. 1972, ApJ 176, 1 
Hauser M.G., Peebles P.J.E. 1973, ApJ 185, 757 
Henry J. P., Arnaud K.A. 1991, ApJ 372, 410 
Holtzman J. A. 1989, ApJS 71, 1 

Hugues J. P., Gorenstein P., Fabricant D. 1988, ApJ 329, 82 
Hughes J. P. 1992, ApJ 337, 21 
Kaiser N. 1983, ApJ 273, L17 
Kaiser N. 1986, MNRAS 222, 323 

de Lapparent V., Geller M.J., Huchra J. P. 1986, ApJ 302, LI 

Lemaitre G. 1933, Ann. Soc. Sci. Bruxelles A53, 51 

Lilje P.B. 1990, ApJ 351, 1 

Melott A.L. 1985, ApJ 289, 2 

Melott A.L. 1987, MNRAS 228, 1001 

Oukbir J., Blanchard A. 1992, AfeA 262, L21 

Peacock J. A., Heavens A.F. 1985, MNRAS 217, 805 

Peebles P.J.E. 1980, The large-scale structure of the universe, 

Princeton Univ. Press 
Peebles P.J.E. 1982, ApJ 258, 415 
Press W., Schechter P. 1974, ApJ 187, 425 
Smoot G.F., Bennett C.L., Kogut A., et al. 1992, ApJ 396, LI 
Tuluie R., Matzner R.A., Anninos P. 1994, Sissa Preprint, 

astro-ph/9403040 
White S.D.M., Frenk C.S., Davis M. 1983, ApJ 274, LI 
White S.D.M., Davis M., Frenk C.S. 1984, MNRAS 209, 27P 
Zeng N., White S.D.M. 1991, ApJ 374, 1 



This article was processed by the author using Springer- Verlag 
LT E X AfeA style file version 3. 



